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ABSTRACT 


Optical scintillation on folded paths has been investi- 
gated for the case of spherical waves and a plane mirror as 
a folding target. A criterion for the statistical indepen- 
dence of the two parts of the folded path was formulated, 
utilizing the path integral technique and the Von Karman 
spectrum. Employing the same approach, an expression 
meet ne scintillation strength (the log-intensity vari- 
ance) on folded paths has been developed. It has been 
shown that the scintillation strength on folded paths 
can be expressed aS a product of the one-way scintilla- 
tion strength and a function of the distance source-detector 
in units of the outer scale. The analytical prediction shows 
that the ratio of the log intensity variances on the folded 


17/6 


paths to the one on the one-way path, goes from 2 £OR 


exact folding (the distance source-detector is zero), to 
77° when the detector is far away from the source. The 


theoretical predictions have been corroborated by preliminary 


field experiments. 
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Pe LN ROD EC TIEON 


Intensity fluctuations of visible or near-visible sources, 
while viewed through the atmosphere (scintillation), play a 
Major role in a large variety of situations such as star 
observations, optical communication channels, Lidars, power 
transmission, etc. It has been understood for a very long 
time that the stars’ twinkling is caused by the inhomogeneity 
of the atmospheric index of refraction. It is only in the 
last thirty years, and especially since the first lasers were 
built about twenty vears ago, that this intuitive understand- 
ing has been put into a rigorous quantitative model. Some 
of the better known contributors to this progress are 
Meearski [1,2}, Fried [3], Lutomirski and Yura [4], Ishimaru 
me, Lee and Harp [6], Clifford [7], Fante [8,9] and Dashen 
[10]. As a result of their work and many others it is now 
believed that in most of the situations the processes involved 
are fairly well understood and the theoretical predictions 
are reasonably close to experimental findings. 

Among the problems which are less understood and for which 
existing theory either does not give reasonable predictions 
or has not been confirmed experimentally, is the problem of 
scintillation on folded paths. In general, this descriptive 
name stands for a class of situations in which EM fields 
propagate in the atmosphere from a source to a target and 


back to a receiver (detector) in the vicinity of or at the 





same location of the source. Typical situations which belong 
menchis class are optical radars, communication links (ground- 
satellite/mirror-ground), remote monitoring of atmospheric 
pollutants, and others. 

The class of point-ahead problems such as pointer-tracker 
(with fast moving targets), power transmission from space to 
earth (which follow a beacon signal from earth to space) 
shares a major question with the previous class, that is to 
say, the correlation between the two parts of the path. 

In this paper we are concerned with scintillation caused 
by propagation along a folded path when the target 1s a 
specular reflector (either a plane mirror or retroreflector). 
Transition cases in which the two parts of the path are par- 
tially correlated are included. To the best of our knowledge, 
enly a smal] number of works in the open literature address 
this problem. We will compare the results of these works 
with ours DUNC reCratead parts Of OU paper. 

In the second part of this’ paper (discussion) we describe 
some accepted formulas, physical pictures and approximations, 
which describe the atmospheric turbulence and the propagation 
of EM fields in it. In the third part we develop an expression 
meme tne correlation between the two parts of the path. In 
the fourth we give the expressions for the scintillation 
strength (the normalized variance of the log-intensity) for 
fapeireine Mirror and retroreflector targets. The fifth part 
describes the experimental part of our work (system descrip- 


tion and results), and a summary constitutes the sixth part. 
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Il. BACKGROUND 


A. THE SCALAR EQUATION OF PROPAGATION 

The effect of the atmosphere on the EM fields propagation 
Manifests itself through the Maxwell equations. The atmos- 
phere is considered to have no free charges (hence zero con- 
ductivity and unit magnetic permeability). It is also 
assumed, based on experimental findings [8], that the time 
scale of changes of the atmospheric index of refraction is 
much longer than the time scale of the EM field fluctuations 
in the visible and near visible band. This assumption allows 
us to suppress the time dependent part of the EM field in 
Maxwell's equations. Based on the above conditions, we find 
that the vector equation for the EM field propagation is 


V*E(r) + k*n*(r)E(r) = V(V-E(r)) ia4) 


where: 


k = w/c is the EM field wave number 


2 
meg GE) Se (5) 
eae thesar electric constant 
Me =) amnoex of refraction. 


It was found analytically [11,12] and confirmed in experi- 


Memes |13), that the depolarization of EM fields in the 


ma 





forward direction due to atmospheric turbulence (variations 
Peis) ) is negligible: Thus the term, on the right hand 
side of Eq. (11-1) which represents the depolarization 
(through mixed derivatives) can be neglected [8], and we get 


the scalar equation of propagation 


v-E(r) + k’n*(r)E(r) = 0 (II-2) 


Bee OLMOSPHERIC TURBULENCE 

As we can see from Eq. (II-2), the effect of the atmos- 
phere on the EM field propagation is caused by the fluctua- 
tions of the index of refraction. The index of refraction is 
a very complicated function of the atmosphere's constitutents 
(aerosols, water vapour , etc.), the absorption coefficients 
[14] (continuous and discrete), the temperature and the pres- 
Sure. Assuming the absorption is negligible (the wavelength 
aoes not match any absorption line or band), the humidity 
gradient is negligible, and the aerosols' contribution to 
forward scattering is small, then the index of refraction 
turns out to be a relatively simple real function of the 


temperature and the pressure [15] 


er een ne) fa + 287245 * ar -3) 
= icy) = \2 
iene : 
P = total pressure in millibars 


iz 





T = temperature in degrees Kelvin 


= 
ii 


wavelength in micrometers 


Even though the expression for the index of refraction is 
a simple one, it is impossible to get an exact solution for 
the scalar equation because we cannot specify the values of 
Eee temperature and pressure along the path. The only way 
we can solve this equation is through the use of their sta- 
mestical properties. 

This approach leads us to express the index of refraction 


in the following form 


meee) = + ny (x) (II-4) 


where No Tomeiceavenoaesvalue Of n(r) along the path and 


n, (x) ftiomiG@eecusSpatial fluctuation. on, is found by replacing 


0 
0! Ty respectively, and it 


obviously cannot cause scintillation. In almost all works 


fee), T(r) in Eq. (II-4) by P 


about the effect of the atmospheric turbulence No 1s taken 
6 Jee ie 

Since the pressure fluctuations (in the microscale) are 
relatively small [16], the fluctuations of n, (x) are governed 
by temperature fluctuations. 

Due to their small order of magnitude (less than one 
degree), the main processes through which temperature inhomo- 


geneities in the atmosphere dissipate, are convection on a 


large scale and molecular dissipation on a very short scale. 


i 





Therefore, temperature inhomogeneities (fluctuations in 
space) which have scale size larger than a certain critical 
value, are closely related to the velocity field within the 
atmosphere. 

Using this relationship between these fields and the 
Kolmogorov model of the velocity field, Tatarski [1] gave a 
Beata stical description of the index of refraction fluctua- 
tions in the atmosphere. 

The Kolmogorov [17,18] model is based on the following 
assumptions: 

ie rhe fluid is incompressible. 

2. The fluid is in a steady state. 

3. The fluid is getting a constant power per unit mass 
from an outside source. 

4. The Reynolds number of the flow (onset by the outside 
source) is much greater than the critical Reynolds 
number of the air. 

5. Kinetic energy dissipation into heat can be done 
only on avery short length scale. 

Since the Reynolds number of the original laminar flow 
is larger than the critical value for air, once the flow is 
Smseurbed, it will turn into a turbulent flow with vortices 
whose length scale is limited only by the dimensions of the 
flow. This scale is called the outer scale of the turbulence 
(Ly), and in the surface layer it is usually taken to be the 


height of the point above the ground. 
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If the Reynolds number of the flow within this vortex is 
larger than the critical Reynolds number, the vortex will 
break down into smaller vortices. Since the Reynolds number 
is proportional to the scale size of the flow, the formation 
of smaller vortices reduces the local Reynolds number within 
each vortex until it is smaller than the critical value, and 
the flow is stabilized. The smallest vortex size is called 
the "inner scale of the turbulence" (0) . 

In the literature about EM fields in the atmosphere, 
these vortices are sometimes called eddies. 

See GNOEMBLE AVERAGING, STRUCTURE FUNCTION AND TURBULENCE 

Speco RUM 

The averaging process through which we get the moments 
of the EM field should be done over many points. Each solu- 
tion of the scalar equation represents a different ny (x) 
field (the ensemble average). Since the n, (x) field has a 
relatively short correlation length, averaging over the 
eeeseeribution of n, (x) along a propagation path can replace 
the ensemble average. In the literature the ensemble aver- 
age is denoted by < >. We will use the same notation. 

Assuming the n, (x) field is isotropic and stationary, 
we can define the following functions. 

1. The autocovariance (sometimes called correlation 


Or covariance) function 


By (4,-Eo) = <n, (4) n,(r5)> (TI=>) 


ike: 





Py. The structure function 
Z 
D (Ey 15) = <{n, (4, )-ny (z,)] > (II-6) 


3. The spectrum function o,, (K) 


\ 


ie LK(X)-E> 
BL {£)°E>) = JJ Je, we ° dk 


Since we assumed isotropic turbulence, the last 


equation can be written as 


B(p) = 4n f $,(K) San K* dk (ae) 
0 


In the last four equations: 


ro spatial vectors 

p a Rap > Bs! 

i —eeSoditdt Frequency vector (of the n, (x) 
field) 

K = [KI 


Using the Kolmogorov model, the spectrum function is 


given by 


= N73 


» i poeseak for 2n/Ly < K < 2m/L, (11-8) 


pt 


16 





Thies ember val isaealled the inertial range. C 

represents the strength of the turbulence and it is 

proportional to the variance of n, ata (OS erice lite 

In most early works, this spectrum was used (by way of 

Peeeoximation) in the range 0 < K < ~. This procedure as 
pointed out by Yura [19], can cause serious deviation in the 
outcome of some calculations. Thus a somewhat different 
spectrum is now being used, namely the modified Von Karman 


spectrum. 


2 ~(K20)* 
-033C, e 


o, (R) a (i= 9) 
ero 


Hill and Clifford [20] showed that a high K values, 
the o, (R) is not a decreasing function, but rather it in- 
creases to a local maximum and then goes to zero as the 


modified Von Karman spectrum. 


mee APPROXIMATE FORMS OF THE SCALAR EQUATION 

The most commonly used approximations for the scalar 
equation are the Rytov (small or smooth perturbation) approxi- 
mation, and the parabolic equation approximation. 

The Rytov approximation is based on the perturbation 
theory and has one basic assumption: 
(IIT-10) 


\VinE Sevan 


1 . 


ay, 





This restriction is equivalent to the assumption that the beam 
can be scattered by turbulence eddies only once, which means 
Bmee it restricts the range of applicability of this approxi- 
mation to weak turbulence levels. 

Starting with the scalar wave equation, and replacing n 


py Ll+n and E by exp(T), we get 


iL 


Z i - 
1 


Vor ee (Vp ees a = 0 


; _ 2 
Taking T = To + Ty and ny — 2ny, we get 


2 2 2 2 2 2 2 
V Ty + (VT) +k° +/Y¥V T, + (VT,) + 2(VT)) (VTL) + 2k ny = 9 


We can see that condition (II-10) is equivalent to 


| VT, | << IVT, | 


a 
ema that if e oe EO is the vacuum solution of the scalar 


equation, the last equation can be written as 


2 2 _ 
asec ia + 2k-n, = 0 (II-11) 


0 is known), for the moments of T,- 


Experimentally it was found [7] that the Rytov approxima- 


moach can be solved (VT 


tion is valid as long as 


af = 124 «7/8 pt1/6 2 go (IT-12) 


n 
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where: 


k = wave number, 
L = distance of propagation, and 
ee = the turbulence strength. 


The Markov approximation is based on the assumption that 


(for an EM field which propagates along the x axis) 


5 du 
5 | eq || Bathe En (II-13) 
OX 
where 
- —~ E(x,y,2) Z 
u Liles, Z)) exp (LKx! (II-14) 


As Tatarski showed [2], the condition (II-13) 1s equivalent 
memene condition l << Kg which is almost always met in the 
atmosphere (Lo ~ 1lmm) for the visible band. 

Inserting (11-14) in the scaiar equation and using condi- 
tion (II-13) gives the following equation, which is called 
the parabolic equation 


0 2 2 rs 
Pik aya (x, Xp) ~ "x, u(x,X)) + 2k n, (x, 0) ulx, > = Cee 15) 


a 


where Xo is a vector in the y,z plane. 
The physical meaning of condition (II-13) is that the 


Bem@ecering due to the eddies, is primarily in the forward 


io 





Peeection (8). Eq. (11-15) can be solved by using the Markov 
approximation [2] or by path integral. Dashen, who was the 
first to show that the parabolic equation could be solved 

by the path integral technique, argues [10] that the con- 
straints imposed on the solution of the parabolic equation by 
both techniques are the same, and identical to the one imposed 


by the parabolic equation. 


E. PATH INTEGRALS AND THE SOLUTION OF THE PARABOLIC EQUATION 
In the forties, the path integral technique was developed 
R.F. Feynman in an effort to use the principle of least 
action as the basic principle of Quantum Mechanics. A good 
description of the technique and its application in Quantum 
Mechanics and a few other fields of physics can be found in a 
feecbOOk by Feynman and Hibbs [21]. The use of the technique 
in solving the parabolic equation is described ina report 
by D.L. Palmer [22]. 
The use of the path integral technique in solving differ- 
ential equations can be demonstrated by the following example 
(free particle in a one dimension motion). 


Applying the Shrodinger equation to this case gives us 


dW n° 52 
i i re (+ boundary conditions) (Ii-16) 
ec 2m 5x2 


where ~ is the probability amplitude function. 
We know that for a free particle in a one dimension 


motion, 
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S = | enmre (une 7) 


Starting from Feynman's basic assumption, namely that 
the probability amplitude assigned to the ith path (connecting 


moints a,b) is 


b 
Vv; (arb) J exp[is(a,b)/n] dz, 


Geeeneche integral is being carried out along the ith path. 
Since the total probability amplitude is the sum of the 


contributions from all possible paths, it is given by 


w(a,b) = fD (paths) exp[is(a,b)/n] (II-18) 


Comparing Eqs. (II-16) and (II-18), we can see that the 
solution of the differential equation (II-16) can be found 
by the use of the path integral technique, providing we can 
come up with the kernel of the differential equation, which 
in our example was obvious. 

In the case of the parabolic equation, the procedure 1s 


more elaborate, and since this is not the main topic of our 


2 





work, we would rather merely quote the result in Section 


Meer .c) of this work. 


F. SCINTILLATION AND SATURATION 

Viewing a source of constant output through the atmos- 
phere, we see that the intensity varies with time (scintil- 
lates), just as in the case of stars twinkling. 

Usually the scintillation strength is measured by the 
variance of the normalized intensity 


2 nl > 


Di Dit - - : 
pee < (1715) > <I/I)> = 1 (Dae) 


where: 


Since it was found that in weak turbulence, the log ampli- 
tude of the EM field is a random variable with normal distri- 
bution [2], the log of the normalized intensity is distributed 
normally in this region. Hence it is useful to work with 


the variable 2% = in(I/I)) and to define the turbulence strength 


by 
a = a 25. ai). (ja 20) 


Defining AI = I-I), we can see that if AI << Ig, since 


<AI> = 0, then 


Ze 





= 
lI 


In(I/I)) = In{(A4t+1I,)/I5] 


ln (1 +S) Sart 
; 0 


Thus Eq. (II-20) can be written as 





2 
ae) ee > Sapp = 4h (II-21) 
: 0 0 2 
0 


There are several models by which the scintillation can 
be explained. The simplest model which was originally used 
by Tatarski [1], describes the eddies in the atmosphere as 
Spherical lenses whose index of refraction differs from their 


Surroundings by An =n These lenses focus or defocus 


1° 
(depending on the ny values) the rays onto the target plane. 
It was believed that this picture gives good approxima- 
tion only in the weak turbulence region (not more than one 
scattering event), until Yura [23] showed that it can be 
applied to the strong turbulence region also, if we incor- 
porate in this picture the notion of the transverse coherence 
length (XQ) and its effect on the focusing ability of the 
eddies. Due to focusing and defocusing of the beam while 
propagating through the atmosphere, the transverse coherence 
length of it is gradually reduced and therefore, the ability 
of the eddies to focus the beam further, is reduced too. 
Hence the growth of the scintillation strength is saturated 


free tt Should reach a maximum. 
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Other models which are used, either in detailed calcula- 
tions or in qualitative explanations, are the phase screen 
model [6] and the scattering model [24]. The first describes 
the effect of the atmosphere on the beam as a distortion of 
the phase front due to different optical paths in the turbu- 
lence. In turn, the distorted phase front will produce an 
interference pattern in Lune target plane. The second views 
the EM fields at the target plane as the sum of fields which 
were scattered at small angles by off-axis eddies and the direct 
(unscattered) field. 

Since in all models, the scintillation strength depends 
on the turbulence strength (density and effectiveness of the 
eddies) and the total number of focusing or scattering events 
along the propagation path, it is found very frequently in 
the literature that the integrated turbulence strength (the 
effective turbulence strength) is given by the Rytov param- 
eter which was defined in Eq. (II-12). Sometimes it is 
denoted by 8 or On: 

It was found experimentally [7] that in the weak turbu- 
lence region (also called the Rytov region), in which the 
mm@eensity distribution is log normal, af goes up to 0.3. As 
of grows beyond this value, the intensity distribution is 
close to the K distribution [25], the scintillation strength 
reaches a maximum [26] and starts to decrease. As of grows 
further, the intensity distribution is best described by the 
exponential distribution [10,27] and normalized intensity 


variance approaches l. 
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tLe Sol tio rhliCAL DEPENDENCE OF THE TWO 
PARTS OF THE FOLDED PATH 
A. PREVIOUS WORKS 

AS we mentioned in the first chapter of this work, the 
question of whether the two parts of the folded path are 
statistically independent, is a basic one when we try to 
describe the statistical behaviour of the EM field which 
propagates along it. Once we can assume statistical indepen- 
dence, we can treat the two parts separately, applying the 
well known results for the one-way propagation to each of 
them. Since our main object is the scintillation along a 
folded path, with a special emphasis on short distances be- 
tween the source and the detector, the first subject of our 
analysis is to try to define a criterion for the statistical 
dependence between the two parts. 

Among the works that we were able to find in the litera- 
ture, which analyzethe folded path question, some like the 
theoretical work by Yura [28] and the experimental work by 
Ben-Simon et al., [29], assume statistical independence of 


the two parts. Yura gives the condition for independence as 


where d is the distance source-detector and Dos Ds are the 
aperture diameters of the source and the detector respectively. 


Since in Yura's paper the target is diffusive (volume of 


aS 





air molecules) and the reflected beam is totally incoherent, 
we cannot apply this result to our case (specular target). 

in other works [30,31,32] the question of the indepen- 
dence of the two parts is related to the wavefront tilt or 
angle of arrival (AOA) fluctuations which depends on the large 
scale eddies. Since we are concerned with scintillation 
which is caused primarily by the small scale eddies, we can- 
not use their results. 

Before we continue our analysis, we would like to mention 
without going into detailed explanation, that the usual cri- 
terion for independence between two diverging paths (their 
end points should be separated by a few isoplanatic patches), 
cannot be applied in our case, due to the fact that ina folded 
path case certain parts of it (near the target) might be very 
close and certainly correlated to some degree. As we will 
see later, this correlation affects the scintillation and 


we Cannot ignore it. 


B. THE EXPRESSION FOR THE PHASE VARIANCE 

We start our analysis with the effect of the atmospheric 
turbulence on the phase of the propagating EM field. 

The phase fluctuations (S) relative to the vacuum value, 
due to the non-uniform medium, are given by the integral 


eeeeng the optical path. 


L 
S = k ” n, (x) dx Cori) 
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Since omar (x)> a) 


L L 
<S> = <k f nj(x)dx> = k f <n,(x)>dx = 0 
0 0 


The variance of the phase fluctuation is 


Tg = <S° >> - <§> a <Sa 
and 
5 L L 
Soe- = <(k J n,(x)dx)(k f n,(x')dx')> 
1 i 
0 0 
5 L L 
ee eee) ice endian Sone (x). (sc?) > 
i e 
0 0 
5 D L L 
soe = ik i eee f es Bee} (IITI-2) 
n 
0 0 
Since in the Kolmogorov model, B(x rx') = BL ty) where 
y = |x-x'| and it is the same function regardless of the 


wergan of y, aces can be approximated by 


5 iy Ae 
oa — wie i BL (y) dy (III-3) 


ai. 





provided L >> y where y is the range at which 
B_(y) << B (0). 

As we can see from Eq. (III-2), the variance of the phase 
is related directly to the field of the index of refraction. 
This is the reason we chose to utilize the phase variance in 
the criterion for statistical dependence/independence between 
the two parts of the folded path. 

In this part of the work we chose to work with an infinite 
plane mirror as a folding target for two reasons: 

1. It is easy to compare the case of a folded path (with 
a plane mirror) when the angle between the two parts 
Sietie Pacem 2s 25, to the case of a straight path 
(28 = 7) (see Fig. 1). There is no discontinuity of 
the rays and there are no singular points at the target, 
as in the case of the corner cube. 

2. When an infinite, or at least very large mirror is used, 
we do not have to take into account the possibility 
that part of the beam will not be reflected. Thus at 
this stage, we can neglect the effect of beam wander. 

melas RELATIONSHIPS BETWEEN THE MEANS OF THE EM FIELD IN 

VACUUM, THE EM FIELD IN A NON-UNIFORM MEDIUM, AND VALUES 

OF THE PHASE VARIANCES 

Following Dashen [10], we can write an expression for ¢ 

eugeee & = Re(ec) at a point (L,r) due to a point source at 


(0,r.). This analytical signal is the solution of the para- 


Ey 
bolic equation [22] and it is given by 


28 








as 
€(0,r),L,r) = 3 {D(paths) 
ee 
x exp{ik f [$(r'(x))“-n, (rte .x)]dx} (III-4) 
0 Loe ere eX 
where: 
dr (x) 
tee) ae 
2 = Suineewector alonguxsascs.. 


From Eq. (III-4) we can get the expression for the analytical 


Signal of the vacuum field by setting n, = O. Thus we get 


0,c,,L,r) = {D(paths) 
Mae ye 2 
ER ff Ile ea) Il gers) (III-5) 
0 


Taking the ensemble average of Eq. (III-4) we get 


—— > 


i 
<5; /D(paths) > (III-6) 


ik (> ee fm 
x exp {= f [x’ (x) ] “dx - ik ,! n, (x(x) +e, x] dx} 
0 


i es ee 
sz /D(paths) x expi> 4! [x' (x) ] “dx} 


L 
x <exp{ik f ny (x(x)+e,x)dx}> 
0 


Zo 





Using Eq. (III-l1) we can write 


L 
<exp{ik f n, (x(x) +e,x) dx}> = <exp(S;)> 
0 
where Si = is? and sF is the phase fluctuation along the path 


P. Since it is generally accepted [33,34] that the phase 
fluctuation in the atmosphere is a Gaussian random variable 
with zero mean, ee is also Gaussian with zero mean, and 


i 
* 
therefore it satisfies the following formula | ) 


eo < (gt) *> 
alt 7 Ji 
a exp [——4——] 


Thus 


* 
( Por Gaussian random variable a, with zero mean and with 
standard deviation o 








ere) co 2 
<e*> = J ef (4,0) da = : i etexp[- + .] da 
—< (2716 -© eG 
20 Cee 2 co Dee 
a = f exp [222 Saya = epic i2) { ie Sn) ae 
Y210 —~ 20 (2710 00 20 


© 


= exp[o*/2] f £(a,0°) da 


=m OO 


exp [o*/2] 


mere £(a0,x) 1S the normal distribution function for mean = x. 
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is” 


P, 2 
<e Dee — exp [- SS io 


Using Eq. (III-2) we get 


L 
<exp{ik / n 


f 1 (E(x) +e,,x) dxt> (III-7) 
2 L L 
= exp{-S ff BULL (x(x)-x(x') “+ (x-x!") 177? ]axdx")} 
0 0 


Assuming that the optical paths in our case, which are the 


main contributors to the path integral, satisfy 
Pace) << |x Sa" | (III-8) 


Eq. (III-7) can be approximated by 


L 
<exp{ik f ny, (xr(x)+e,x) dx}> (III-9) 
0 
Ke L L 
= exp{- a i f B(x, x!) dxdx'} 
0 0 


M@eerting Eqs. (111-5), (I1I-9), (IfTI-2) in Eq. (III-6), we 





Set 
0 <S7> 
se(0,r,,L,r)> = € (0,r,,L,r) exp[- —5—] (i =e) 


Syl 





where <s*> is the variance of the phase fluctuations, due to 


propagation along a path of length L. 

From Eq. (III-10) we see that the effect of the turbulent 
atmosphere on the average of a propagating EM field, can be 
expressed by the variance of the phase fluctuations along the 
geometrical path. 

It is obvious that in vacuum, apart from a phase shift 


due to the reflection, the fields 23 and Ear, as described in 


Fig. 1, are the same. We can ask ourselves what will be the 


relationship between <En? and <E5,>? Since in the case of 


meame> there is no folding at all, it is obvious that in this 


2L 


case there is no statistical dependence (relating to the 


folding). Hence, the closer Sen? Ee SEnr7 7 the less depen- 


dence (statistical correlation between the two parts), we 
have. 
Using Eq. (III-10) and the observation that up to a phase 


- = er: we come to the conclusion that the statis- 


tical dependence between the two parts of the folded path, 


can be determined by the relationship between <s<> and <S5_>: 


the closer their values are, the less the two parts of the 


Semstant ¢€ : 


folded path are correlated. Therefore, in the next steps we 
will calculate 35 and i (we omitted the angular brackets 
for the sake of convenience), compare their values, and formu- 
late a criterion for statistical independence. 

Before we go to the next step of the analysis, we have 


to evaluate the effect of the constraint which was formulated 


in Eq. (IIfI-8), on the applicability of our approach. 


BZ 








Figure 1. The reflected (folded path) and transmitted 
EM fields at distance 2L from a point source. 
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PemGomenempelmlr (%,(x)) 15 a point On the optical path 
through the turbulent medium (which 1s given by the normals 


to the wavefronts), the constraint (III-8) is equivalent to 


dr (x) 
dx 


which means that the normals to the wavefronts are almost 
parallel to the x axis, which in turn means that only small 
angles of scattering are allowed. This constraint is the 
same as the one for the utilization of the parabolic equa- 
tion [8,10] (Eq. (I1I-13)), and both are met for almost all 


Semceivable conditions. 


D. THE VALUE OF i. IN TERMS OF s*, FOR THE VON KARMAN 
SPECTRUM 
From Eq. (III-2) we can write 
2 a Ls 
a ee ae Be > x" )) 
0 0 
and 
> > _2L Dt 
ee CK a ax ) dx' B (|x-x!'|) 
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w= xe xX 5 See ot 
See il 
JF) = aD: 
Li, 2L.—-Z 2 tne Z 
3 «2 at aL 0 a 
o> ats «odz J B(jzjjdy+ f daz f B,({z|) dy} 
i 0 Z i —Z 
7 
i 0 
E el dz B_(|z|) (2L,-2z) + - Jf dz BA ([z]) (2b, +22) 
i 


where L. is either L or 2L. 


Taking z > -Z in the second term, we get 


ie 
2 2 7 
So = 2k° ff (Lb, - 2)B, (z)dz 
a 0 
Inserting the explicit forms of BA fs) {2] and the Von Karman 
spectrum for 9, (K) 
z . Sin (Kx) 42 
B(x) = 47 ) OG cca) mere K~ dK 
-033 Cc“ (x) 
SL es) renee ene 
n (K* + 15%) 77/8 
Assuming C* (x) = constant and taking y = KL, we get 


35 





L 
;/h 


20 0 (Li/Lo72) ay y sin(zy) 
ee) az 5 f oe UG 
a 0 0 (y +1) 
where: 
A = .033 22> 1 k* ce 1272 
n “0 
and Zap ey is 


0° 


Being formula 2.3(11) of Ref. [35] we get 


L,/Lp 
2 :, ~ 7 3 
Sy. (AY i; dz(L,/L, Zz)Z Ky 3 {2) 

a 0 

where: 
eos ee (Com Piao) 
m, 70 
and K (Zz) = Modified Bessel function of the third 
1/3 ao 


Solving the integral in Eq. (III-ll1) (detailed solution 


is given in Appendix A), gives for (L;/Lp) 2 


s aya /*r (5/6) 
"L. = ; eee 0 (ITI-12) 


which is the same as the expression which is given by Tatarski 
[2]. Thus we can see that for large (L/L). 3 is approxi- 
a 


mately linear in L. and 
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S = 2S (Lit 13) 


Z 
i 


BE. THE EXPRESSION FOR S 
The geometry of the folded path and the notations which 
we use in this section are described in Figure 2. We assume 
that the mirror is large enough so that: 
1. Diffraction effects can be neglected. 
2. There are no fluctuations in the amount of energy being 
intercepted by it. 
We also assume that the field at the point D (E,) can be 
found by the path integral technique, using SMD as the geo- 


metrical optics path. It should be pointed out that since 


the calculation of the Ee value by the path integral assumes 


oe << ] (scattering at small angles) 


the application of the technique to this type of path might 
be questionable, for the reason that it includes one point (M) 
at which the above condition is certainly not satisfied. 
Our justification for the use of the path integral tech- 
nique along this type of path, is based upon two arguments; 
1. In the case of fields propagating in vacuum, we know 
that apart from a phase constant, ED = Eos: Therefore, 
we can use the path integral technique (Eq. (III-5)) 


along the path SMD by excluding the point M from the 


integration. 


om 





Figure 2. 





O 


(y,X,) 
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= 
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The geometry and notations of the folded 
path case. M = a plane mirror; S =a Dolne 
Source; D = detector. The distances SM and 
DM are of length L. 
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2. If the results of our analysis in the two extreme 


cases; 
a. 9 = 0 
be c= sa 4 
match the results which are obtained in other ways, 
the fact that at the point M, or tx) is infinite, has 
no effect on the applicability of the path integral 


technique to our case. 


At the end of this section we will discuss the second 


argument. 


As we saw in Section III.C, in order to find Ze we have 


to calculate sé by 


st = x? ff axdx' BL(|x-x'|) 


4 SMD 


where the integration is along the line SMD. Thus 


L D vane 
J dxdx' B ({x-x'|) 
L n 


Kk? J dxdx' B (|x-x'|) +k 
0 n 


te] RD 


22 





where integration on the interval 0,L means integration along 
SM, and on the interval L,2L means integration along MD. 


Pielmceqrals IL. and I x' is the distance of the point 


3 4! 
(y,X,) from the point M along the line MD, and the notation 
|x-x'| stands for the actual distance between the points, 


which is denoted by P. 


Using the transformation y = x-L; y' = x'-L 
9 1s; 
yk ,) f aydy' B(ly-y'|]) = 2, 
Since BL(|x-x']) is symmetric in x, x', ee and 
therefore 
2 ha 
me = k leepegciecaiba | Be (il xe= 2x2 ||) - (III-14) 
0 n 
5 L 2L 
in d bray =x = 
s x 2 dx' B (|[x-x']) I, +1, 


eon Eq. (III-2) 


mmenevaliuation of I, 


From Fig. 2 we can see that 





x' = (L-x,) /Cos (28) ; adx' = - dx, /Cos (25) 


(III-16) 
y = (L-x, )Tan(26) 
pe x" | = p = [ (L-x, ) “Tan? (28) + (x, -x) 2177? (III-17) 
Hence I, can be written as 
Ke L L 
I, = cos (3B) ) dx f dx, B (P(x, x,)) 
2 Sin® 


Applying the transformation (X, xX, ) + (P,a) and observing that 


ee, = 1; = P Sin(29+a) 





Sini(Ze ) 
2. xX, = L- P Sina Cot{(28) 
ox 9X 
X)Xy dP 00, 
3 Jls—5 ) = ae 550 = P Cot(2¢d) 
—it _t 
dP da 
meeror 0 < a < 71/2 - 8 
0 < P Leon (2c) P, 
Sin (28+a) 
eeror 71/2-8 < a < 1-280 
QO < p L Sin(2o) = Pp 
= ad Sale ae. 2 
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we get 


acepee eee ()) dP 


+ f da f PB, (P) dP} 
1T/2~-8 0 


Making the transformation 8 = 1-(28+a) in the second 


integral, we can see that it is equal to the first one. Thus 


2 n/2-6 SL 
eS ere da f PB_(P)dpP 
3 Sin(2é) n 
0 0 
where 
> - L Sin(29) 
ae Sin(28+a) 


Inserting the explicit forms of mee) and the Von Karman 


(RK) and assuming that o, (Kr) = o,(K), we get 


1/2-8 yy 00 a 
me «(CA f da f Shy Gk 2(z*+1) 11/63 in (yz) 
0 0 0 
where: 
Wieceenc Cc. 
oS n 0 
Sin(28) 
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fees anc ©} 


Y. L Sin (26+a) 


Using formula 2.3(11) of Ref. [35], we get 


1/2-6 ¥1 
= 4/3 
I, = A, J da J Gye ys? 
0 0 
where: 
Ors3 22/3 2 ke ee oe 
na Sin(20) T(ll/6) 
K = Modified Bessel function of the third kind 


al 
it 


the gamma function 


Comparing the last equation with Eq. (III-12) gives 


22/3 I, ; se" is 4/3 
oe UL a Ae ey) 
: Sin(26)n-/*r(s/6e)L = o 0 1/3 


Inserting the last equation and Eq. (III-15) into Eq. 


(III-14) gives 
s“ = 28°[1 + £ly.,8@)] 
1 oY 


iS 


Sieoy UuSing (III-13) 
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2 
eee cae fy, 8) 1 (a= 108) 


where 
y 
A 1/2-86 at 
ee ae 4/3 = 
ey?) = 5 i da f y Ky /3(y) dy (III-19) 
0.0 0 
and 
7 - 52/3 
2 od Ona 7G) 
7 fap een ae) 
0 Ly 
ne Yo/Sin (26 +a) 


New sOlved Eq. (IIi-19) numerically and found out that for 
small angles 9, which is the case in most of the applications, 
the £ dependence on 6 is very weak, and therefore 
Elyy,8) = fly). 

The result of the numerical integration is described in 


Figure 3, and we can see that 


lim Elyg-9) = 1 (when 6 > 0, viaual 0) 
ama’ 
0 
lim Elyg-9) = 0 
G>1/4 


and therefore 
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9 (6 Trews 203) °K Zo 
uotqzouny Be Se (6T-III °5a) (9’-A)F FO ANOTAeYSq UL 


Bt 
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8 
S 
> 
€ 
3 


NOTLONAA ADNSGNAdad SHiL'd 


°¢ oANbTty 


a 


(X) + 


é° 
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2Sor 69 = Q 
S<(L, 8) a (III-20) 
Z Es 
Sor, 6 = 7/4 


These are the values of = for the two extreme points of 


i 
mae transition region (from exact folding to statistical 
independence of the two parts of the folded path). In 
previous works on the subject of scintillation on folded 
paths, the authors assumed either exact folding or statistical 
independence of the two parts of the folded paths. Thus, as 


we pointed out in Section III-E, we can compare our results to 


those calculated by others, only at these two points. 


ieee CRITERION FOR STATISTICAL INDEPENDENCE 

We have now enough information about the behaviour of 
feet, 3), in order to define a practical criterion for statis- 
tical independence of the two parts of the folded path. 
Strictly speaking, the two parts are independent only for 
9 = 7/4. Yet because Elyg,8) decreases very fast at small 
Yo and has a very long tail (approaches slowly to zero), we 
can define the region at which we have statistical indepen- 
dence, as the region in which L(Y 979) < ae From Figure 3 


we see that this criterion is met (for small 98). 


ae 2S 
Since 
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Meece sin (28) 
Y — ae 
0 Lo 


for small angles 6, Yop can be written as 
eae“ o 


where Dis the distance source-detector. Thus we see that the 
statistical independence between the two parts of the folded 
path depends on the distance source-detector, expressed in 
units of the outer scale. 

As we said at the beginning of this section, one of the 
justifications for employing the path integral technique in 
our case, even though the main constraint on its utilization 
1S violated at the point M, is the observation that the re- 
sults in the extreme cases (98 = 0; 6 = 71/4) match the results 
which are obtained by other techniques. 


For the case of @ = 71/4 we got 


2 


fin = 7 
Sea) = Se, 


which is the expected result based on physical arguments and 


the results of Section III.D. 


For the case of 9 = 0 we got 
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When a beam propagates back and forth along the same turbu- 
lent medium, the phase fluctuation due to each path element 
Ax is doubled. Thus comparing the random variable S in the 
case of a path of length L (S,) to the case of a folded path 


with 8 = 0 (S.), we get 


ry 
tC bw 


which fits our results. 
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Pee oer N PeEATIONSS CRENGTH FOR THE FOLDED PATH 
ee LHe RATIO BETWEEN THE SCINTILLATION STRENGTHS FOR THE 
FOLDED PATH AND THE ONE WAY PATH 
The log intensity variance for the Rytov region for 


spherical waves is given by [1] 


1/6 JENS: 


2 = ie 
To (Ep = a ebek oe L (IV-1) 
In the strong turbulence region (exponential distribution 
of I), the scintillation strength is given by the variance 
of the normalized intensity [8] 
ad - 69 
Tr = Ll + — aos 275 (IV-2) 
(o,) 


Comparing Eqs. (I1II-12) and (IV-1) we get for the Rytov 


region 
2 bo DOr. 2 yo z 
a, (L) Rye G yo L SF A(k,L))L Sy (Iv-3) 
TT k Lo 
or for a path of length 2L 
2 z ByeOr 2797 One a 
Cy (2h) = A 2 L Sor (Iv-4) 


Prem Eq. (IIlI-18) we get for the folded path 


NO 


Po(r) = A 5ii/é 3/6 


Z 
L S, li + Ply qr) ] 


~ 
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We can express the ratio between the folded path scintilla- 


Z 
iv 


(97 (L)), in the weak turbulence region, by 


tion strength (0, (F)) and the one way scintillation strength 


oO” (F) 
31/6 





(ee Elyy/9) J (IvV-5) 


NO 


o, (L) 


This ratio is the parameter which we are going to measure in 
the experimental part, in order to verify our model. The 


dependence of Ron Yo (for small 8) is described in Figure 4. 


B. DISCUSSION AND OBSERVATIONS 


ie) one Range OE Reana the Saturation of Scintillation 
Sienengien 


As we can see from Eq. (IV-5) and the corresponding 


figure, the range of Ris 


51/6 oi7/6 


Therefore when we work with a system for which 9 is close to 
zero or zero, we should expect the scintillation strength to 
be about seven times greater than the scintillation strength 
of the one-way path. 

One should be aware of this large ratio especially 
while designing a system or an experiment which is based on 
the log-normal distribution. Since the region of the log- 


normal distribution, as we mentioned in Chapter II, is given 


by 
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faethe case of exact folding (6 = 0), the upper limit of 
this region can be reached along relatively short distances 
and/or relatively weak turbulence levels. 

On the other hand, since saturation can be achieved easily 
(relative to the one way case), it provides a way of carrying 
out experiments in the saturated region without the difficul- 
ties of using very long paths. 

2. Comparison with Other Works 

MihGmelso, 57) tOlund that the scintillation strength 
of a spherical wave which propagates along a folded path 


(plane mirror) with 6 = 0 is 


Zee 2 K°x (2L-x) 


1G; oO 
2 2 , 
G{F) = 641° k - dx 4 OC ote ak | (EV=G) 


Since he used the phase screen approach [6] and the Kolmogorov 
spectrum, his results are applicable only in the Rytov region 
and there is no dependence on the outer scale. 

Using the Kolmogorov spectrum and the Mellin trans- 


mem (formula 6.5(15) of Ref. [35], we get from Eq. (IV-6) 


L 
2 _ _ <i 2i= x) = 5/6 
rece Cl) a= AS ! | 


a2 





where 


176 Ceo. 32 


A = .033(4)2 iaayeyt—1) Cos(sn/l2)n «k Co 


3 


Making the transformation 


y = xX/L 
we get 
2 iy a 5/6 
were 2 = 0) = A, L ee yiviG2—y) 1 (IV-7) 
0 
It is well known that for the one way path of length 
il 2 | 
o2(n) = 16 12k? C ax f Ko (K) Si 2K (=x) | ax 
5 ; n = 2kL 


Following the same steps as in the previous case, we get 


L 
2 _ x =x) / 6 
ie a 


where 


A, = .033 4(-1) [(-5/6) Cos(51/12) qe es - 


Sie) 





and we can see that 


A, = 2/76 A, 
Taking y = x/L we get 
b2 Ve oy 13/6 ; 
oe es, mi salva enya (IV-8) 
From (IV-7) and (IV-8) we get 
of (Fé = 0) a3 J t(2-y1/6ay 


a 
a, f ty(l-y)1?/Say 
0 


2 
O, (L) 

Taking y = 2x in the upper integral and observing that the 
integrand in the lower integral is a symmetric function about 


y = 1/2 we get 


ila 


277P a 2 4/6 f [x (t=) 197 Pax 


R = ee = 5177/6 
2 5/6 
A, 2 f [x (1-x)]~% “ax 
0 
which is exactly the same results as we got by using the path 
integral technique. 
It is easily seen that even though the results for 


8 = 0 are the same (using the two techniques), Smith's paper 


cannot give the values of R for 8 # 0. 
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Gamo et al., [38] measured the value of R (for 6 = Q) 
along a path of 274m at night (using a collimated beam HeNe 
laser), to be 6.668. From their paper, it is not clear 
whether the measurements of of (L) and o* (F) were done 
simultaneously. 

TweeonceWetghne Funectron an the Case of Exact Folding (0 = 0) 

The relative contribution to the scintillation of 
different points along the path is described by the weight 
momet1On. 

From Eqs. (IV-7) and (IV-8) we can see that the weight 
functions W, (x) for spherical waves are: 


For one way Wo(x) = [x(l-x)]?/° 


For the folded path Wi(x) = eee) oe 


The two functions are shown in Figure 5. 

From Figure 5 we can see that in the case of the folded 
path (@ = 0), the main contribution to the scintillation comes 
from the vicinity of the mirror. In the case of the single 
path, the main contribution is from the center of the path. 
Since it is very unlikely that we can find the proper condi- 
tions under which C* (x) along an atmospheric path is constant, 


we should be careful when we compare results from the two types 


Of path. 
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4. The Effect of Path Folding on Aperture Averaging 


The discussion and the evaluation of the scintillation 
strength in previous sections assumed that we have a point 
detector. The term point detector refers to a detector with 
aperture diameter (D 5) smaller than a critical length (Dp) - 

The aperture averaging can be explained qualitatively 
in the following way [7,40]. Due to atmospheric turbulence, 
the EM field at the detector plane has a certain spatial struc- 
ture. This implies that the intensity (the modulus of the EM 
field) has a spatial structure too. In order to measure the 
full effect of the scintillation utilizing direct detection 
of intensity, the detector aperture should be smaller than 
the distance along which the intensity variations are negligi- 
ble (Dy), otherwise the detector will average different 
intensities across its aperture. (In the language of Young 
[40], the detector will filter out the high spatial recente 
cies of the intensity's spatial variations.) The critical 
length is given by 


_ ~1 
B,(D)) = B (0) e 


where BL 1s the covariance function of the intensity at the 


detector plane. 


It is well known [1,2] that in very weak turbulence 


(az) 1/2 


Dy ~ J = 
where \ is wavelength and L is distance of propagation. 


a 





In the case of strong turbulence, it is not clear 


what is the exact value of D, [8]. The difficulty in obtain- 


0 


ing the D, value arises from the difficulties in calculating 


0 
mme fourth moment of the EM field. Since the transverse 


correlation length of the EM field (r [23] is a decreasing 


0? 


function of the distance whilst den 1S an increasing one, at 


a certain distance the field and the intensity spatial covari- 


ance functions are dominated by r Thus, beyond certain 


0° 
Euebulence strength, Do should be a function of Io: 

As Yura [28] pointed out, the MCF of the EM field 
1S proportional to the phase variance. 

It follows from the above mentioned arguments that, 
as the phase variance increases the critical length Dy 
decreases. Since we found that in the case of the folded 
path 


2 


ee = ee) + aia) ise 


it is obvious that while working with the folded path, we 
should utilize detector apertures which are much smaller 


than those which are used in the one way case. 
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V. EXPERIMENTAL WORK 


The main object of the experimental part is to measure 
the values of the function R(y 9) and to compare them with the 
analytical function which is given in Eq. (IV-5). 

In order to obtain the values of Rly), we have to 
measure simultaneously the values of the one-way scintillation 
strength (9%) and the scintillation strength for the folded 
path (0%). These measurements should be done at several source- 
detector distances (assuming that L, is the same for all 


0 


measurements). 


Ewe ios EXPERIMENTAL SYSTEM 
The experimental system consists of two major parts: 
1) The optical system; 
2) Data reduction system. 
mae ene Optical System 
The optical system is described in Figures 6 and 7. 
The main problems which we dealt with during the 
planning and construction phases of our work were: 
a) Selecting proper sources. 
b) Making sure that the beams which propagate along the 
one-way path and the folded path see the same turbulence. 
c) Obtaining high quality optical elements with rela- 
tively large apertures (within our budget and time 


Span constraints). 
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The sources of our optical system should meet certain 
requirements: 

a) Low noise. They should have low values of intensity 
variance in order not to obscure the intensity fluc- 
tuations due to the atmospheric turbulence along 
relatively short paths. 

b) Operation ina pulsed mode or a simple way of 
modulation. (This requirement evolves from the 
way we were planning to subtract the background 
noise.) 

c¢) When passing through the same turbulence, the respec- 
tive measured scintillation strengths should be highly 
correlated. 

d) Synchronized operation either in pulsed or modulated 
modes of operation. This requirement makes it possible 
to compare pairs of pulses. 

Technically it is easy to employ a GaAs laser at the 
plane mirror end (the one-way source) in a pulsed mode, and 
a mechanically chopped HeNe laser at the other end (folded 
path source). These sources in these modes of operation, 
were used before by the Electro-Optics group of the Physics 
Department at the NPGS in previous experiments [41]. 

dee lne =Sounces 

It was found (assuming the sources' intensities 

are normally distributed) that the normal intensity variances 


of the sources were: 
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3) 


HeNe laser (Spectra Physics Model 155) mechanically 


~4 


chopped; ae oe 0 


i 


12 elements array of GaAs; S : 


aon 


8x8 elements of GaAs; sé oo oa This laser 
matrix which is manufactured by Laser Diode Labs, was 
equipped with a Lucite parallelopiped rod mixing ele- 


ment, in order to enhance its output uniformity. fThis 


mixer was designed and built at the NPGS. 


The conclusions from these measurements were: 


i ) 


i) 


In view of the expected levels of scintillation strength, 
all three lasers can be incorporated as sources in the 
optical system. 

The 8x8 GaAs 1S a better (more stable) source than the 

12 elements array GaAs. Therefore we decided to continue 
the experiments with the 8x8 GaAs laser only. ‘In the 
following text, the term "GaAs laser" refers to this 
laser. 


The correlation coefficient between the HeNe and 


the GaAs lasers were measured in a configuration which is 


described in Figure 8. The GaAs in a pulsed mode was triggered 


by the mechanical chopper via a pulse generator. In this way 


we were able to get interlaced series of pulses with a con- 


trolled time separation between them. For each laser we 


measured the log-intensity variance (in a way which will be 


described in Chapter V.B) and each pair of variances (one 


for each laser) was compared. For each run we calculated the 
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statistical correlation coefficient between the two series 

of variances' values. One run of this part of the experiment 
is described in Figure 9. In this particular run we got the 
highest value of correlation coefficient (.973). In the other 
runs the correlation coefficient values were much smaller. 

As a result of this measurement, we decided to 
replace the GaAs laser with a second CW HeNe laser (similar 
to the first one), which is modulated by a Pockels cell 
modulator. (We tried to use a second mechanical chopper in- 
stead of the Pockel cell modulator but failed, because we were 
unable to synchronize it with the first one.) We found that 
the values of the two series of scintillation strengths (log- 
intensity variances) of the two HeNe lasers were highly 
correlated. A typical run with o9 = .996 is described in 
Figure 10. 

The measurements of the correlation coefficients 
were carried out in the basement corridor of Spanagel Hall 
at the NPGS. (Its structure and turbulent features are des- 
cribed in Ref. 41.) The path length was 145 meters and the 
time separation between the pulses of the two lasers was 
eoo=—350 us. 

Based on the results of these experiments, we 
decided to continue our experiments with the two modulated 
CW HeNe lasers. 

b. Frozen Atmosphere 
As we mentioned in Chapter II of this paper, the 


atmosphere can be considered frozen for time intervals of 


or 
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Peeco lb ms..oince the reund trip time of flight along the 
folded path (when the target is up to 10 km from the source) 
is less than 70 us, the atmosphere can be taken to be frozen. 

In our experiments we employ two interlaced 
series of pulses with a time separation (At) within each 
pair of pulses. The time separation in our system is dic- 
tated by the equipment which was available to us, and is not 
less than 250 us. 

We tested the assumption of frozen atmosphere 
by measuring the correlation coefficients of two series of 
log-intensity variances. The two series of pulses originated 
from one source by uSing two series of slits on the mechani- 
cal chopper for the HeNe laser, or by the use of the double 
output signal of the pulse generator, for the GaAs laser. 

For Ehewceuble pulse” HeNe laser experiment 
with t = 270 us, the correlation coefficient was op= .998. 
This experiment which was carried out ine the corridos” shows 
that for time separation of up to 270 us, the atmosphere is 
frozen. 

This experiment served also as final checkup for 
the equivalency of the two channels of the data reduction 
system (see Chapter V.B). 

Cc. Optical Elements 

We Originally planned on doing the experiments 
with two targets; a plane mirror and a corner cube. The 
Original configuration of the optical system was dictated 


by four major requirements: 
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1. High quality elements. To maintain the phase relation 
between different parts of the beam. 

2. Large aperture elements. To intercept a large part 
of the beam, in order to avoid the effect of different 
transmitted/reflected intensities on the scintillation 
Siegen til. 

3. The laser at the target end (the one-way transmitter) 
cannot be placed behind the target, because its beam 
cannot be transmitted through the corner cube. 

4. The beam-splitters in use should not cause an inter- 
ference pattern. 

The Graginal configuration of the optical system 
is described in Figure ll. 

The beam-splitters and the plane mirror target 
are made out of large pieces of flat fused Quartz, which were 
coated with Al, to achieve the required reflectivities, at 
the NPGS facilities. They were checked for flatness by 
interferometric techniques. A full description of the testing 
and coating processes is given in Reference 42. 

The interference patterns due to the double reflec- 
tion from the two surfaces of the beam-splitters, can be 
eliminated by applying an anti-reflection coating to one of 
the surfaces. Since we were not able to use this technique 
for various reasons, we used a different approach. 

The beam-splitter at the detector end is located 
very close to the detector (5-10 cm). Hence, by using thick 


plates (taking apart the two secondary sources) we can eliminate 
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this problem when the beam has a small angle of divergence. 
This beam-splitter thickness is 1/4", it has an elliptical 
shape (4 x 6"), and was coated to give a ratio of 1:2. The 
beam-splitter at the target end is a circular (6" diameter) 
plate of 1/2" thickness. AS we can see from Figure ll, its 
interference pattern (at the detector) was eliminated by 
applying the Al coating to the back surface and working with 
a polarizer and angle of incidence equal to the Brewster 
angle, to eliminate the reflection from the front surface. 
In the course of the experiments, it was found 
by observing the beams at the detector end, at night, that the 


beam which originates at S, 1s very homogeneous but the 


2 
homogeneity of the beam that originates at Sy suffers a great 
degradation due to its double pass through the Quartz beam- 
splitter at the target end. Sinceweuse a back surface 
reflection for the first beam (it passes this beam-splitter 
twice too), the only explanation for the different effects on 
the two beams, 1s that the inhomogeneities in the beam-~ 
splitter material are of larger length scale than the diameter 
of the beam near the laser output aperture. 

Because of the beam quality degradation, we changed 
the optical configuration from the one which is described in 
Figure ll, to the one which is described in Figure 6, uti- 
lizing the circular beam-splitter which has 70% reflectivity 


as the target mirror (with the coated surface facing the 


detector end). 


as 





As a result of the above mentioned change, we 
abandoned our original plan to measure the R values with a 
retroreflector as a target. 
d. Other Elements 
1) The mechanical chopper is constructed out of circular 
metal plate which is rotated at about 3000 rpm by an 
electrical motor. There are 20 slits on the circum- 
ference of the metal plate with a ratio of 1:4 between 
the slit's width and the opaque regions between them. 
The chopper modulates the CW HeNe laser beam at a fre- 
quency of about 1 kHz with "pulse" width (laser on) of 
about 200 us. An LED (GaAs) and a detector which are 
placed on opposite sides of the chopper plate, provide 
a trigger signal at the same frequency as the laser 
modulation for the whole system. Due to this arrange- 
ment, we do not have to use a stabilized chopper. 
2) The Pockels cell modulator is made by Coherent Asso- 
Clates (Model 3003). 
3) The diameter of the aperture in front of the detector 
is 3 mm. 
Ze Data Reduction System 
The data reduction system is described in Figure 12. 
It 1s a modification of a system which was used by the EO 
group for measuring the scintillation of one source. This 
System and its mode of operation are described in detail in 


Bets. {41, 42]. 


We 
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The elements which were added by us are: 

1. The second channel (demodulator and log converter). 
2. The electronic switch. 

ee 6CC LOCK. 

4. Software (for the Nic-80). 

The electronic switch (AD 7510D1) alternately sends 
the signals from the parallel channels to the computer. 

The clock is constructed of a set of pulse generators 
which are triggered by the mechanical chopper's LED-detector 
pair. The time sequence of the signals is described in 
migure 13. 

The change in the software enables us to store the 
measured log-intensity values of the two series in different 
storing spaces, for later calculation of the log-intensity 
Variances. 

The main steps which are carried out in each channel 
of the data reduction system, are the following: 

1. The Demodulator samples the signal intensity and later 
(between pulses) it samples the background intensity. 
It subtracts the second from the first and sends this 
difference to the log-converter. 

2. The log-converter takes the log of the input signal and 
sends it to the computer input via the electronic switch. 

3. The computer samples the input signal, converts it from 
analog to digital and adds one to the proper counter. 


After sampling 16,384 pulses, it calculates the variance 
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assuming that the distribution (after the conversion 

to log intensity was made by the log convertor) is 

normal. The calculated value is printed on the 

assigned printer. 

Note: When we measured the sources' (lasers') normal 
intensity variances, we bypassed the log con- 
vertors, and used the average value of the 
intensity distribution to normalize it. 

After we completed the modifications, we tested the 
two parallel channels for equality by comparing two series of 
pulses which came from the same laser. The high correlation 
coefficient between the variances of the two series (9 = .998) 
proves that the two channels are equal. The results of this 


experiment are described in Figure 14. 


B. EXPERIMENTAL RESULTS AND DISCUSSION 

The final experiment was carried out along a path of 62 
meters on the lawns of the NPGS. The height of the optical 
path above the ground was 1.1 meters. 

Because of weather conditions, we made measurements on 
February 27th and on March 3rd and 4th, 1982. Throughout the 
experiment, the weather was fair to clear. There were light 
winds (0-3 knots) from various directions, large portions of 
the path were shaded by high trees, and most of the time the 
grass was wet. The meteorological parameters which were taken 
from the NPGS meteorological station (measured at the NPGS) 


with other relevant parameters are summarized in Table I. 
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Information about the conditions under 
which the final experiment was done 


Distance Source R Day Time Tr RH% Weather 
Detector (cm) 


0 Ga 3-4-82 1237-1245 56 63 Clear 
20 6.19 2-27-82 1626-1633 58 80 Scattered 
clouds 
40 Gee 2 3-3-82 1447-1454 56 67 Scattered 
clouds 
70 sles: 2-27-82 1239-1247 61 68 Clear 
O'S Boe 2-27-82 1053-1100 58 69 Scattered 
elouds 


During all the runs, the wind was between 0-3 knots from 


wemrous directions. 
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Figure 15 shows the experimental results, and the function 
R(Yg19) Weenie 1s Given in Eq. (IV-5). 

The R value of each experimental point is an average of 
16 measurements. Each measurement is the ratio of the log- 
intensity variances of the folded path (given distance source- 
detector) to the one way. 

The Yo values of the experimental points were calculated 


by 


oP eso 
where D, is the distance source-detector (for the folded path), 
and Ly was taken as 20 cm, to given an optimal fit with the 
analytical expression. 

Even though we have only five experimental points, it 
seems that the experimental results match the analytical 
Meearclion quite well and confirm it. 

The outer scale value (Ly ~ 20 cm) which gives an optimal 
fit between the experimental results and the analytical pre- 
diction, seems to be relatively low, for a path which is l.l 


meters above the ground in an unstable surface layer. For 


a stable surface layer, this is a reasonable value [44]. 
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VI. SUMMARY 


Optical scintillation along folded paths with different 
angles between the two parts of the path, is a part of the 
noise in many systems which measure the intensity of a re- 
flected signal (e.g., a system which measures pollutant 
concentrations in the atmosphere via their absorption). 

In previous works the scintillation strength along these 
types of paths was calculated by assuming either that the 
two parts of the path are statistically independent or by 
working with an "exact folding” (9 = 0). 

In this work we utilized the path integral technigue and 
the Von Karman spectrum, to formulate a criterion for the 
statistical independence of the two parts (for spherical 
waves and a plane mirror). This criterion was given in terms 
of the distance source-detector expressed in the outer scale 


unit, namely 
vine = D/Ly > 4 


For the same case, we developed an analytical expression for 
the ratio (R) of the log-intensity variances (folded path 


to one way path) as a function of Yor in the Rytov region 


3i1/6 


Rly) - hie Ely,)] 


where 
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0 Ye = 2L(L = the distance source-target) 
f(y.) =d¢et y, > 4 
0 0 — 
A ijumax QO (exact folding) 


and showed that this prediction gives for the end points 

(Yo = 0+Y5 = 2L) the same results which were predicted in 
other papers which dealt only with these points. Our ap- 
proach gives an analytical prediction for intermediate values 
of Yo: which others do not. 

In order to verify our predictions, we designed and put 
together an experimental system by modifying an existing data- 
reduction system, and constructing the optical system. This 
Beem cnabled us tO measure simultaneously the scintillation 
strength along a one-way and a folded path, which to the 
best of our knowledge, was never done before. 

We measured Ryo) ana tound that 

1) the maximum and the minimum values of Rly) fit the 
analytical prediction well. 

2) the slope of the experimental value fits well the 
predicted one, for an outer scale value a0 =e Ole 

Based on these results, we think that our analytical 


expression predicts well the scintillation strength on folded 


mechs. 
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APPENDIX A 


THE SOLUTION OF EQUATION (III-11) 


maeenight—-hand side of Eq. (III-11) 


L,/L, 


ile, 3 
S =a f dz(L,/L)-2)z 4 Ky 2 2) 


is a difference of two integrals. The first one is 


er 
if sere ues ft ms az zt/ sx (Zz) 
a. Pa 0 0 eS 
and the second one is 
m= A Pik Zen) 
2 a 73 


Baewsirst integral 
Using formula 7.14(6) of Ref. [43], we get 


~2/3 1/2 


ci = A, (L,/L)) 2 [ZH (A-1) 


1 ato 672 UK 


Ws -2/3 ‘2) 
L/h, 
ee 


Hy 73 '2)B_9 38211 9 


where H(z) is the modified Struve function. 
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noucerts 7.2 12) 7.2013) and 7.5(57) of Ref. [43] give 


the following relationships. 


wo 


I (2) = ) (2/2) eM*E vat p (mer+l) (A-2) 
= m= 0 
TT 
K (2) = Sanaa ar - I, (z)] (A-3) 


H(z) = J (2/2) *@*** op (me3/2tr) TP (m+3/2) 17} (A-4) 
m= 0 


For z > 0, formulae (A-2) and (A-4) can be approximated by 


(2/2) ~/T (r+1) (A-5) 


tH 
N 
il 


i (2/2) tt tp (3/2) 7 (3/2tr) 17+ (A-6) 


N 
il 


Inserting (A-3), (A-5) and (A-6) in Eq. (A-1), and calcu- 


lating the lower limit gives 
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= 1/3t/? lim{az-b gon - ae - Ae 
Z>0 
where a, b, c and d are constants. 
Thus for z = Q 
21K) ,,(2)H_> 73 (2) = Hy 73 (2) Ky 3 (2) | = 0 (A-7) 


Deemetatrge zZ, formulae 7.13(7), 7.5(52), 7.5(63) and 7.13(4) 


of Ref. [{43] give, up to mae 
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T(r - 5) 


(A-8) 
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H(z) = SE py - EGA) (92) 71 + 0 (277) 1 (A-9) 
- V27 r(>+ 4) 
il: 
N-1l T (m+ =) 
— =f J (-2)™ ——* (2/2) 572 F 02 ANI) | 
m=) [(r+5-m) 


Using Eqs. (A-8) and (A-9) in the upper limit of Eq. (A-1l), 


we get (z = L,/L, Soa) 


(2)H_4 7362) + (z)K (Zz) ] 
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— nye {1 + Flee} te! (L;/L,) 
+ O((L./L.)~*)} Gea 
ee) 


Pienwsecond integral 
Using formula 7.14(3) of Ref. [43], we get 
Li/lLy 


4/3 


i = A, [-z Ry 3 '2)] 


0 


For the upper limit we make use of Eq. (A~-8) and we get, up 


Zoe) 


-2 
ZO Z (L/L, 


4/3 
Z 3 2) 


(By 1229/8 e721) + eed (22) + 0(z*)] 


mame for L./Ly = 2 2 10 is almost zero (less than OMT 


For the lower limit, by use of Eqs. (A-2) and (A-5) 
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fein es, , (2) = sma Men ey saae 
Z>0 Z +0 
TT gl/3 


Sin(4n73)F (173) 


meas Cor L/L, > 10 


A, T 2i/3 
as cea ore ly on (A-12) 


We can now write the solution of Eq. (III-11) 


ey 2 
(L,/Ly)1 Tr (5/6) 


sh A a RRR - Hr ytg 
+ O(L,/L)) *) - a 
Sr, _ Bel Sp  agitg) i (FIPS Z a1) ones 
> i/2 en 0((L,/Lp) 2) } 


~ Sin(4a73) 1 (17311 (576) Mi Uo! 
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Ae aie dt x8 


NUMERICAL EVALUATION OF Ely)78) 


The equation for a) is 


A 1/2-6 ¥y 
- 4 4/3 
f(y,-2) = Joo da y) y°’ “Ky /3(y) dy 
As we showed in Appendix A, 
M4 y 
ik 1 
4/3 4/3 
fey Ka.tyidy = <-y “~K, ,4(y) (1) 
0 173 4/3 0 
i738 
me 4/73 _ 1 2 
7, Y "Kay3'Y) = sin(aa73)P(-173) 2 
Thus 
A T™/2~9 
2 _ 4/3 
oye) = 7, J da [B-y)" “K, ,3(y,)] 
where 
3 
T 2 
- Sin(4n73) 1 (-173) oad 
Or 
ip 
£(y,,8) = oe - Telyy,9) ] (B=2) 
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where 


1t/2-6 


Tg(¥ge) =f da yy Ky /3(¥y) 

and 

Y, = Y9/Sin (28 ta) 
Taking 

8 = 26+a 
we get 
I. (y)/8) - a AByy/ 7Ky /3(¥}) . st a8 yt! "Ky 73 (¥,) 
From 

Y; = Y)/Sin 8 
we get 

Yo Cos 8 
heist aoe of 
Sin 5 
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Therefore, for the first integral in I 

















5 
ie sin*g| 9%) 
Cos 8 Yo 
and for the second integral 
Fn Sing oval 
Ges s Yo 
Since 
2 
ae Sin’3 _ +0 
Cos 8B (1-sin?g) 1/? , Gg eye 
1 eat 0 
Meye,?) = y J en en | nay, 
a 0” 0 | aa 2) 1/2 0 ’ oan 2, 1/2 
where 
oe Sin(26)/L) (B-3-a) 
ee = Yo)/Sin(26) = L/Ly (B-3-b) 
= - 2h Sah aae, 
Y = y,/Sin(t/2-8) = y,/Cos§ = ——=— (B-3-c) 
b 0 0 Ly 


Integrating by parts (making use of formula 7.14(3) of Ref. 


[43] gives 


oe: 





vew..4/ 3 ve 
oy K (y) - i 
ae dy = See Livsy yp Oe (y) 
Zz Daye y 0 4/3 
y/ Ma = VG) 0 Yo 
Q 
Y 4 
+ f Sec™*(y/y)) y°/ 7K, ,3(¥)] 
YQ 
Since 
Seq ~ (1) = Q, 
Sec l(y Vee 1/7 2-29 
ay 0 y 
sec ty yi) = § 
bi-<0 : 
and 
4/3 
y R43 'Y) 
has a finite value, 
—— : Ay Ss 4/3 
I. (Yq) = Cin 2 BOY Ky 3 Yq) rae Ky 73 'Yp) 
Ya Yb 
Sic) Cy) Cy alg Gy) dy 
YQ YQ 
where 
g(y) = sec” (y/yy)y°/ "Ky »3(y) - (B=4) 


o2 





Hence 


A 
E(y),6) = = [B (1 /2-8) + (n/2-0) ya’ 7K 


(y,) (B-5) 


4/3 


y y 
ec f aiyidy = | gyal 
Yn Ky /3(¥p , g(y) dy , gly 
0 0 


This the exact form of Elyy,0). 


From the definition of y (Ec. (B-3-b)) we can see that 
a 


4 


in all practical situations, when the path of propagation is 
parallel to the ground but not too high, we can take ae i 


Beang tormulba 7.13(7) of Ref. [43], we get for y >> l 


4 —m aS 
y er, ,,(y) : ss Oty Ne 


pay AE ETD LAME Br oy 


Pe eee ess 7a) Olly 2) 


(iy 2) 


By taking only the first two terms in the brackets, the error 


4/3 


in the value of y Ky 3 (Y) Pomeess tiiatieels for y > 10: | Ehus 


temey > 10 we can take 


4/3 


ye’ Ky gly) 2 (a/ eae SS 72 


Dene ; 


(an) 1/2y5/6 o-¥ 


JA 


Comparing this expression with B, we find 


a8 





a3 


B/1000 > y Ky 3) Eor 10 


tS 
|v 


, 4/3 
Simce in Ba. (B-5), B and ee iy ey have the same 


= 


coefficient (1/2-8), we can approximate this equation by 


A 
Elygr8) = =AIB(n/2-8) + ey,’ 
0 


Ya Vy 
- f gly)dy - f  gly)dy] 
Yo Yo 


Numerical Evaluation 
The numerical evaluation of Eq. (B-6) was done in the 
Computer Center of the NPGS, on an IBM 3033. The program is 
written in FORTRAN IV and is given at the end of this 
appendix. 
The important information about the program is listed 
below. 
ie the function E(y9,8) is evaluated with Yo as the 
independent variable. The step by which Yo 1s changed, 
1s determined by a constant (parameter) of the program. 


2) The angle § is calculated from the equation 
5 
9 = 5 Sin (L./Y) 


3) The parameter L, which is given for each run, is the 


il 


distance source-mirror in units of the outer scale. 


Lu; = L/Ly 
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4) The constants Ao, Beacesconstanes Of Ehe program. Their 
values, which were calculated from Eqs. (III-19) and 


(B-l), are 
A = ./7934 Be =) el. 250 


5) The values of the modified Bessel function of the third 
kind K(y), are calculated in different ways for the 


two regions: 


Gymmore << 4 


(ob) 4 < y 
iMmecegLom (a) we used formula 7.12(21 of Ref. [43]. 


oO 


ac) { Cosh(rt)exp[-y Cosh(t) ]dt 
0 


which can be written 


=O 


av = f  cosh(rt)exp[-y Cosh(t) ]dt + R(r,t,) 


Using 
Cosh(t) = =(e +e") > > e 


aD 





we get 


Rules en) ee - J (eee ne tC) exp(- a) at 
0 —- 2 2 
1c 
0 
Taking 
Es 
x = aa BG = e e = dt and Se seit 
0 x 0 
we get 
R(r,Xp) << 5 i (xt 44 (Fol) Wem? a5. 
X 
0 
Boo tr = 4/3> re=l = 1/3 and -(rtl) = -7/3. 
Therefore 
R(4/3,xX)) < = J pooja! cas 
Xx 
0 
“X,y/2 rns = 
ee 0 2 0 
= me 5 ye iy 
ong 
“Xp y/2 
= 
R(4/3,X9) < - (Xp pel ayy) 
BOmem — 7s? E-l = -2/3 and -(rtl) = -4/3, 


Therefore 
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we =A é 
Riggss,x.) < 1 foe xY/2 3 es” 
0 — 2 
Xo y 
since in bq. (B=6) 


and we calculate E(Y5 19) Onilye LOE 2 2, Seeting 


iG, 4) 107’ gives 


0? 


Xo “ee SUmLOr © = 473, and 
Xo wl Om ret. © = 1/3. 
Thus, we conclude that 
6 
KS) { Cosh(rt)exp[-y Cosh(t) ]dt (B-8) 
0 
+ R(r,6) 
where 
R(r,t) < 10°’ for xr = 4/3 or 1/3, and y > .l. 


We calculate K.(y) ECOmin co eDyPIncegrating Eq. 
(B-8) numerically in steps of At = Mop FOEeYy - a4, 
we use formula 7.13(7) of Ref. [43] 


N 


= il 
K.(h) = (n/2y)*/*e7¥[ b em sSt i (2y) 0 (y§)] (3-9) 
n= 


oi 





We take N = 10, and therefore for y > 4, the order of 


the absolute error is less than 


7 = oon = 
ay ae a On = 10 8 


op For Y, OF Y,z which are greater than 30, we use the 


approximation 
Ya'Yb 30 
i g(y)dy = f  gly)ay (B-10) 
Y avs 
0 *0 
For Y, > 30 we can write 
Yy SHO) 
J ogty)dy = Jf gly)dy + Rly), 
YQ YQ 
where 
a! 1 /3 
in ~ 4 
Riy,;) = 40! Sec “(y/yg)y “ "Ky ,3(y) dy 
"L473 
Soy 2) re y Ky 3 (y) dy 
4/3 “2 
= tn/2) [oy °K, 73) 
0 
4/3 


ete Your oO, | ¥ Ky7/3 ¥) is a decreasing function. 
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Thus e 


R(yy) << (4/2) (y, ~ 30) 30° 7K, ,3 (30). 
2 4 
ifaeengys0 < ye < 10 , 
Geyer) <= 1/2)’ 10 30°7°e 29 =~ 107! 


7) The numerical integration of fg(v)dy is done in steps 
of ay = 201. 
Additional information concerning the numerical 


integration can be learnt from the program itself. 


The Program 
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